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Recent fluorescence spectroscopy experiments on single wall carbon nanotubes reveal substantial 
deviations of observed absorption and emission energies from predictions of noninteracting models 
of the electronic structure. Nonetheless, the data for nearly armchair nanotubes obey a nonlinear 
scaling relation as a function the tube radius R. We show that these effects can be understood in a 
theory of large radius tubes, derived from the theory of two dimensional graphene where the coulomb 
interaction leads to a logarithmic correction to the electronic self energy and marginal Fermi liquid 
behavior. Interactions on length scales larger than the tube circumference lead to strong self energy 
and excitonic effects that compete and nearly cancel so that the observed optical transitions are 
dominated by the graphene self energy effects. 
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The optical transition energies of semiconducting nan- 
otubes, along with their dependence on the nanotube di- 
ameter and chiral angle have been studied in a recent 
series of fluorescence spectroscopy experiments Q, H, |^ . 
Though the experiments were originally interpreted in 
the context of a simple non interacting electron model, it 
has become increasingly clear that electron interactions 
play an important role in determining the optical tran- 
sition energies 0, H, S 0- pointed out in early work 
by Ando|j], interactions lead to (1) an increase in the 
single particle energy gap and (2) binding of electrons 
and holes into excitons. More recently, Spataru et al.0 
have reached a similar conclusion by computing the opti- 
cal spectra for selected small radius nanotubes. However, 
the systematic dependence of the transition energies on 
nanotube radius has not been addressed. 

In this paper we examine the optical excitations of car- 
bon nanotubes in the limit of large radius, R, where they 
inherit their electronic structure from that of an ideal 
sheet of two dimensional (2D) graphene. This permits a 
systematic study of the radius and subband dependence 
of the excitations to leading order in . In this limit 
the electron interactions fall into two categories: (1) ID 
interactions on scales longer than the tube circumference, 
and (2) 2D interactions on scales smaller than the tube 
circumference. We find that the ID long range inter- 
action (1) leads to both a substantial enhancement of 
the energy gap and a large exciton binding energy which 
both scale as Although both effects are large they 

have opposite sign and ultimately lead to a moderate en- 
hancement of the predicted optical transition energy. By 
contrast, we find the 2D interactions (2) lead to a log R/ R 
correction to the bandgap renormalization. This singular 
behavior can be traced to the effect of a the Coulomb in- 
teraction on the dispersion of 2D graphene, which leads 
to marginal Fermi liquid behavior 8]. This logarithmic 
correction is not cancelled by the exciton binding energy, 
and leads to a nonlinear scaling dependence of the tran- 
sition energies on R. The presently available optical data 



indeed show this non linear scaling behavior and agrees 
favorably with the predictions of the large radius theory 
even for tubes with moderately small radii i? ~ 0.5 nm. 

Below we review the non interacting electron predic- 
tions for the energy gaps of semiconducting tubes and 
show that they can not explain the nonlinear scaling be- 
havior present in the observed transition energies. We 
then present the theory for large radius tubes, focusing 
first on the effect of the 2D interaction on scales shorter 
than the circumference. We then incorporate the longer 
range ID interactions into the theory. 

The simplest model of nanotube electronic structure, 
based on non interacting electrons in a linear graphene 
spectrum, predicts that the energy gaps of semiconduct- 
ing nanotubes are 



E!1{R) = 2nhvF/iR, 



(1) 



where R is the nanotube radius, n = 1,2, 4, 5 de- 
scribes the 1st, 2nd, 3rd and 4th subbands, and vp is 
the graphene Fermi velocity. For a tight binding model 
on a honeycomb lattice with with lattice constant a and a 
nearest neighbor hopping amplitude 70, tivp = V37oa/2. 
The linearized model (1) is exact in the limit of large 
radius, and is the first term in an expansion in pow- 
ers of 1/R. Corrections due to curvature Q and trigonal 
warDing|inj are proportional to 1/ sm39/ R^ , where 9 is 
the chiral angle = denotes an armchair wrapping) 
and v = ±1 is the chiral index. A central prediction of 
the non interacting model is thus that for large R the 
band gaps scale linearly with n/R - a, fact that can be 
traced to the linear dispersion of graphene at low en- 
ergies. The large R limit is most accurate for nearly 
armchair nanotubes for which the sin 30 corrections are 
smallest. For such tubes Eq. 1, describes the tight bind- 
ing energy gaps to better than 1% for tubes with radii 
as small as 0.5 nm. The next term in the expansion 
at 0{1/R^) is neghgible. Here we focus exclusively on 
nearly armchair nanotubes, where large R scaling can be 
meaningfully applied, sin 30 corrections, when present in 
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FIG. 1: Optical transition energies in the first two subbands 
for semiconducting nanotubes measured in Ref. 1 as a func- 
tion of n/3R. The filled/open symbols correspond to [p,q] 
nanotubes with chiral index v — p — q mod 3 = +/ — 1. The 
dashed line is prediction of the non interacting theory. The 
solid line is Eq. (4), which incorporates the effect of the 2D 
Coulomb interaction. 



specific nanotubes, lead to deviations from the scaling 
predictions 6, j|. 

The observed transition energies do not obey this hn- 
ear scaling behavior. Extrapolated to large radius, the 
ratio between the first two transition energies appears 
to saturate at a value of 1.7, rather than 2 - a fact we 
have called the "ratio problem" Q. In addition, the ob- 
served transition energies are systematically larger than 
the non interacting prediction. A nearly armchair nan- 
otube of radius 0.5 nm has an observed first subband gap 
of 0.98 eV, which corresponds to Fermi velocity of 7.35 
eVA. This value is significantly larger than the value 5.3 
eVA found in graphite and the value 6.1 eVA deduced 
from resonance Raman data. The "blue shift" cannot be 
represented by a simple scaling of the transition energies 
since it is larger for larger radius tubes and therefore not 
linearly proportional to n/R. 

In Fig. 1 we plot the transition energies reported in 
Ref. 1 as a function of n/3R, where n is the subband 
index, and R is the tube radius deduced in Ref. 1 by 
exploiting the pattern of sin30/i?^ corrections. We have 
used different symbols to represent the data with j/ = ±1. 
The separatrix between the data for the positive and neg- 
ative 1/ locates the data for the nearly armchair tubes 
with 6* ~ 0. At the separatrix the sin 39 /R^ corrections 
are absent, so that the radius dependence should be de- 
scribed by the large R limit to order 0{1/R^) . It is clear, 
however, that even at the separatrix, the linear scaling 
relation (1) is not satisfied. Nonetheless, it is striking 
that the data near the separatrix for the two subbands 
lie approximately on the same nonlinear curve. The sim- 
plest interpretation of this apparent scaling behavior is 



that these energies probe the dispersion of 2D graphene 
at a wavevector qn = n/3R. This suggests that the ratio 
problem and the blue shift problem have the same origin. 

Gonzalez et al. have shown that the Coulomb inter- 
action in 2D graphene leads to a singular correction to 
the electron self energy. Consider the Hamiltonian 

J I 2 J |r-r'| 

where -0 is a Dirac spinor with two copies for the K- 
K' degeneracy, and n — tp^tp. The Coulomb interaction 
is characterized by a dimensionless interaction strength 
g = e'^/hvp- In lowest order perturbation theory the 
electronic dispersion is 

E{q) = hvFq[l + {g/^)log{A/q)], (3) 

where A is an ultraviolet cut off of order the inverse 
lattice constant. The nonlinear behavior as (7 —> is 
a consequence of the long range singularity of the 2D 
Coulomb interaction V{q) = 27re^/g. It is thus impor- 
tant to account for screening. The semi-metallic Dirac 
spectrum of graphene leads to a static polarizatibility 
Il{q) = {l/4:)q/vF- The linear dependence on q exactly 
cancels the 1/q singularity of V{q), leading to a multi- 
plicative renormalization of the interaction analogous to 
screening in a 3D dielectric. The g — s- logarithmic cor- 
rection to E{q) survives screening although its coefficient 
is renormalizcd. In a static screening approximation the 
renormalizcd interaction is gscr = 5/(1 + 9''^ l"^)- 

Though it is derived in for small g, this result has 
deeper implications, since it shows that the weak in- 
teraction limit is perturbatively stable. (3) is invariant 
under the renormalization group (RG) transformation 
A Ae~^, g — > g{t), vp vpii) with 

dg/M^~g'^/A\ dvp/di ^ vpg/A.. (4) 

We may interpret (3) in terms of a scale dependent renor- 
malization of Wi;^ and g. The interaction vertex = Tivpg 
is not renormalizcd, so that the scaling is characterized by 
a single parameter g. Eq. (4) shows that g is marginally 
irrelevant: at long wavelengths g becomes smaller and 
perturbation theory becomes better. This implies that 
even for strong interactions the system flows to the per- 
turbative limit at long wavelengths where (3) and (4) 
are valid. Therefore, the dispersion for small q is given 
exactly by (4) with renormalizcd parameters vp and g, 
which depend on the cutoff scale A. We thus have a sit- 
uation similar to Fermi liquid theory, where low energy 
quasiparticles behave like non interacting particles, al- 
beit with renormalizcd parameters. Here, however, the 
marginal irrelevance of g leads to logarithmic corrections 
which do not disappear at low energies. As emphasized 
by Gonzalez et al.|8|, this singular behavior is a signature 
of a marginal Fermi liquid. 
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Though (3) is exact for 9 ^ 0, it remains to deter- 
mine the values of the renormahzed parameters vp and 
g when the bare interactions are strong. Using the Fermi 
velocity of bulk graphite (where 3D screening eliminates 
the logarithmic sing ularity), hv^^b.?, eVA we estimate 
a bare interaction strength of go — e'^/TiVp — 2.7 at a 
cut off scale Aq of order the inverse lattice constant. A 
crude estimate of the renormahzed parameters may then 
be obtained by extrapolating (4) to strong coupling using 
g = {g~^ + {l/4) log Ao/A)"^ For A ~ .5nm"^ this gives 
g ^ 1.1 and hvp 12.9 eVA. A more accurate theory 
requires knowledge of the form of the RG flow equations 
(4) for strong coupling and requires an approximation. 
Gonzalez et al.jSl have developed a GW approximation, 
which incorporates a dynamically screened Coulomb in- 
teraction. A simpler theory can be developed within 
a statically screened approximation. We find that the 
results ag ree within 5% with the dynamically screened 
theory[ll|. For static screening the renormahzed dis- 
persion has the same form as (3) with g replaced by 
ffscr = ff/(l + 57''/2). The RG flow equations are simi- 
larly modified with a factor of (1 + gTr/2)~^ on the right 
hand side of (4). This leads to a refined estimate of the 
parameters at A ^ .5nm^^: g = 2.0 ; hvp = 7.2 eVA. 
The screened interaction is gscr = 0.48. 

The nonlinear scaling form of the separatrix in Fig. 
1 is consistent with Eq. (3). Choosing the scale 
A = 0.5nm~^, the data is well fit with the parameters 
vp = 7.8 eVA and g = 0.74. These parameters are in 
acceptable agreement with the statically screened theory 
described above, given the theory's simplicity. The 2D 
interactions in graphene appear to explain the nonlinear 
scaling of the data in Fig. 1 and thus resolve both the 
ratio problem and the blue shift problem. 

Nevertheless the agreement between the data and the 
interacting theory of 2D graphene is surprising because 
the latter does not account for excitonic effects, which 
are known to be large 0, 0, 0]- To describe excitons it 
is essential to account for the ID interactions on scales 
larger than R. In addition to exciton binding, these in- 
teractions lead to an increase of the single particle energy 
gap. To address this issue we have numerically calculated 
both the single particle and particle-hole gaps. We find 
that the two ID interaction effects largely cancel one an- 
other, so that the R dependence of the particle-hole gap 
is ultimately well described by the 2D theory. We will 
begin by discussing our numerical calculation. We then 
show how these conclusions can be understood within a 
simple ID model. 

We have computed the the single particle and particle 
hole energy gaps for nanotubes in a statically screened 
Hartree-Fock approximation. Our calculation is similar 
to that previously reported by Ando|^], though here we 
focus on the R dependence of the energy gaps. We use a tt 
electron tight binding model, which includes all ID sub- 
bands. To avoid the complications associated with the 
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FIG. 2: Single particle gaps (dashed line) and particle hole 
gaps (solid line) for the first four subbands of semiconducting 
[p,p] nanotubes with phase shifted boundary conditions cal- 
culated for 5 < p < 25. The thick line is the prediction of the 
2D theory Eq. 4. 

i/s'm39/R^ corrections we study semiconducting tubes 
by calculating excitations of armchair tubes with an ap- 
propriate phase shifted boundary condition impose an 
energy gap. The single particle band gaps are computed 
by evaluating the exchange self energy using a statically 
screened Coulomb interaction. The particle-hole gap is 
determined by numerically diagonalizing the Schrodinger 
equation for the particle and the hole in the renormahzed 
bands bound by the screened interaction. This is equiva- 
lent to solving the Bcthe-Salpeter equation in the static 
screening approximation. 

In Fig. 2 we plot the single particle and particle hole 
gaps as a function of radius and subband index. To em- 
phasize the corrections to linear scaling behavior we pro- 
vide a log-linear plot of En{R)/E^{R) as a function of 
R/n, where E'^ is given by (1), and is proportional to 
n/R. The prediction based on the statically screened 2D 
theory of graphene given in (4) is shown for comparison 
[1^. The single particle gaps are strongly enhanced rel- 
ative to their non interacting values, while the particle- 
hole gaps are only moderately enhanced. Thus, most 
of the enhancement of the single particle band gap is 
cancelled by the electron hole interaction that binds the 
exciton. Moreover, since the slopes of all of the curves 
is the same in Fig. 2, both the single particle and the 
particle hole gaps exhibit the same logarithmic increase 
with radius. The excitonic binding energy, which is the 
difference between the two, does not have the logarith- 
mic increase, and scales inversely with R. The particle- 
hole gaps for the different subbands lie nearly on a single 
straight line, close to the prediction of the 2D interacting 
theory. This is consistent with the scaling behavior in 
the experimental data in Fig. f . In contrast, the single 
particle gaps are well above the predictions of the 2D 
theory, and do not obey scaling with subband index. 
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The essential features in Fig. (2) can be understood 
within a simpler model for the ID interactions on scales 
larger than the tube radius. For example, consider a 
semiconducting nanotube with a bare energy gap 2A 
with an infinite range interaction V{x) = Vq. This is 
the constant interaction model, familiar from the the- 
ory of the Coulomb blockade. In this model the interac- 
tion energy is VoN'^/2, where N is the total number of 
electrons. The single particle energy gap is then simply 
2A-|- Vq. The particle-hole energy gap, which determines 
the energy of optical transitions is 2A. Since the ex- 
citon is electrically neutral, its energy is unaffected by 
the infinite range interaction. For this model the exciton 
binding energy exactly cancels the enhancement of the 
single particle gap. 

Though the ID Coulomb interaction Vo((7) = 
2e^ log qR is not truly infinite range, the infinite range 
limit is an appropriate starting point for describing 
the ID effects. In the static screening approximation, 
^scr(g) = Vo{q)/il + VoiqMq)). Since the ID polar- 
izatibility Il{q) ~ q^R^/vp for small q, the q ^ part 
of the interaction is unscreened and screening suppresses 
only the shorter wavelength components of the potential. 
This leaves a screened interaction which is more strongly 
peaked at low momenta qR << 1, i.e. closer to the in- 
finite range limit. Note that this is consequence of the 
one dimensionality of the nanotube and has no analog 
in a 3D semiconductor, where the long range interaction 
is uniformly reduced by the dielectric constant. These 
considerations help to explain the behavior in Fig. 2. 
The net effect of the long range ID interactions on the 
excited states is relatively small in spite of the fact that 
the renormalization of the single particle energy gaps and 
the binding energy of the electron hole pair are separately 
quite strong. 

The scaling of the exciton binding energy, Eb with 
R may also be considered in a simple ID model. Be- 
gin with the Hamiltonian (2) defined on a cylinder of 
radius R, and integrate out the high energy degrees of 
freedom down to a cutoff scale A ~ 1/R. The renor- 
malized Hamiltonian then has the same form as (2), and 
depends only on three parameters, e^, vp and R. It fol- 
lows that the eigenvalues of H have the scaling form, 
(e^ / R) f(q), where g = e'^/hvp is the interaction at scale 
Since g is scale dependent, the absence of the 
logarithmic correction of Eb in Fig. 2 implies that the 
scaling function / is independent of g. Perebeinos et al. 
have found an approximate scaling relation for the ex- 
citon binding energy for interactions screened by a dielec- 
tric constant 4 < e < 15. For nearly armchair tubes with 
effective mass m ^ 1/{vfR) they find /(g) ^ g^^^/e", 
with a ~ 1.4. This describes the crossover between the 
Wannier limit e 3> 1, where f{g) ~ g/e^ and a strong in- 
teraction limit e 1 where the dependence on g is weak. 



Note that this scaling argument does not imply that 
the the band gap renormalization and exciton binding 
scale like n/R. The apparent scaling behavior for the 
particle hole gaps in Fig. 2 is a consequence of the cancel- 
lation between the ID interaction effects. Because of this 
near cancellation, the effects of the two dimensional elec- 
tronic interactions can be seen clearly in the experimental 
data. It is interesting that theory presented derived from 
the leading order contributions in 1/i? to the excitation 
energies provides a good description of the data over the 
range of experimentally measured tube radii. 

We also note that the large single particle gaps shown 
in Fig. 2 are likely to be important for many nanotube- 
derived devices, but have yet to be measured directly in 
experiments done to date. They are accessible in princi- 
ple by measuring the activation energy for transport in a 
semiconducting tube, or by measuring the threshold for 
photoconductivity following optical excitation into the 
lowest subbands. Interpretation of the gaps measured 
in scanning tunneling spectroscopy are complicated by 
screening effects from the substrate, and make it difficult 
to extract the single particle gap of individual tubes. 
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